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j» In this contribution, I summarise the studies of the properties of Bose-Einstein condensed systems 

composed of up to twelve pions or kaons carried out by the NPLQCD collaboration. These 
investigations have provided precise determination the 1=2 Tin and 1=1 KK scattering lengths 
and the first determination of three-hadron interactions from QCD, finding a repulsive three- 
Q-i pion interaction of size consistent with naive dimensional analysis and a three kaon interaction 

c | consistent with zero. We have also determined the isospin (strangeness) density dependence of 

the isospin (strangeness) chemical potential, finding results in surprisingly good agreement with 
y—t chiral perturbation theory. 
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1. Many body lattice QCD 

Lattice QCD has had major impact in many aspects of particle physics phenomenology and in 
describing the spectra and structure of single hadrons. Computing resources and lattice algorithms 
have reached a stage where it is now worthwhile to consider more complicated hadronic observ- 
ables such as those in the baryon number, B > 1 sector — the realm of nuclear physics. Here there 
are many observables that are phenomenologically important to nuclear structure and interactions 
and to nuclear astrophysics about which very little (or nothing) is known experimentally or theo- 
retically. Systems containing n > 2 mesons are also of interest in a number of areas from RHIC to 
neutron stars. These systems present a significant opportunity for contributions from lattice QCD. 
Recently, the first attempts to study systems of more than two hadrons have been made by the 
NPLQCD collaboration [1, 2, 3]. The results of these studies are summarised herein. 



2. Multi-meson systems 



It has long been known how to exploit the volume dependence of the eigen-energies of two 
hadron systems to extract infinite volume scattering phase shifts [4] provided that the effective 
range of the interaction, r is small compared to the box size L (since r ~ ra~ 1 for most interactions, 
this constraint is m % L S> 1). In recent works, this has been extended to systems involving n > 2 
bosons [5, 6, 7] and n = 3 fermions [8] in the case when the relevant scattering length, a, is also 
small compared to the box size. The resulting shift in energy of n particles of mass M due to their 
interactions is 
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where the parameter a is related to the scattering length, a, and the effective range, r, by 
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The geometric constants that enter into eq. (2.1) are 
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and "C m =n\/m\/(n — m)\. The three-body contribution to the energy-shift given in eq. (2.1) is 
represented by the parameter r\ 3 , which is a combination of the volume-dependent, renormalization 
group invariant quantity, Tjf , and contributions from the two-body scattering length and effective 
range, 

where 

=H3(M) + ^(3V3-4tt) log (ML) - %S<m • 

The quantity T73 (^i) is the coefficient of the three-7T + interaction that appears in the effective Hamil- 
tonian density describing the system [7]. It is renormalization scale, fi, dependent. The quantity 
5? is renormalization scheme dependent and we give its value in the minimal subtraction (MS) 
scheme, J^ms = -185.12506. 

Lattice QCD measurements of these energy shifts allow one to extract the parameters a and 
T] 3 . To determine the energy shifts, we study the correlators (specifying to the multi-pion system) 

C n (t) oc <(X>-(x,o) (?r + (0,0)) ) , (2.4) 
where n + (x,t) = u(x,t)Ysd(x,t). On a lattice of infinite temporal extent, 1 the combination 

G n {t) - j^y n ^ ®f e-*> > , (2.5) 

where AE n is the energy shift appearing in Eq. (2. 1). 

To compute the (n\) 2 Wick contractions in Eq. (2.4), we note that this correlation function can 
be written as 

C n {t) oc ( ( Jfnr) f ) , (2.6) 

where 

n = £ 5(x,f;0,0)5 t (x,f;0,0) , (2.7) 

X 

and S(x,£;0,0) is a light-quark propagator. The object (block) IT is a 12 x 12 (4-spin and 3 color) 
bosonic time-dependent matrix, and r\ a is a twelve component Grassmann variable. Using 

{ - ai -a 2 _ -a,,^^ ^ ^ £ai a 2 ..a^..^„ gftft ^ ^ , (2 . 8 ) 
leads to correlation functions 

Cn{t) = e**"**^ e m . Ml ..^_ n (n)J (n)£ .. (n)* . (2.9) 

'Effects of temporal (anti-)periodicity are discussed in Ref. [3]. 
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For an arbitrary 12 x 12 matrix, A, 

det(l+AA) = JL e<w e m .. ft2 (1+AA)5(1 + AA)A...(1+AA)S 

+ A" 12 C ;1 e**-**-*^ e^..^..^ ( A )g (A )£...( A )J 
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where in the last line we identify the matrix A with n. Further, 



det(l +AA) = exp (Tr [log [ 1 + XA] ] ) = exp Tr 
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Therefore, by equating terms of the same order in the expansion parameter A in Eq. (2.10) and 
Eq. (2.11), one can recover the n-7i + correlation functions in Eq. (2.9). As an example, the con- 
tractions for the 3-7T + system are 



C 3 (f) - tr C)S [n] 3 - 3 tr c , s [n 2 ] tr c ,s P] + 2 tr c , s [n 3 



(2.12) 



where the traces, trc,s> are over color and spin indices. Contractions for n < 12 mesons ai - e given 
explicitly in Ref. [2]. 



3. Two- and three- body interactions 



The NPLQCD collaboration have computed the n pion and kaon correlators in the previous 
section using domain wall fermion [9, 10] propagators on various ensembles of MILC 2+1 flavour 
rooted staggered gauge configurations [11] (parameters are shown in Table 1 and further details 
are given in Refs. [1, 2, 3]). In order to correctly calculate these correlators for large n, very high 
numerical precision is necessary (our calculations use the arprec library [12]). By performing 
a correlated fit to the effective energy differences extracted from these measurements, we have 
determined the two- and three-body interactions. The two body interactions extracted from this 
analysis agree with those extracted from the two-body sector alone [13]. The resulting three body 
interactions are displayed in Fig. 1. The three pion interaction is found to be repulsive with a 
magnitude consistent with the expectation from naive dimensional analysis. In contrast, the three 
K + interaction is consistent with zero within somewhat larger uncertainties. 
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Table 1: The parameters of the MILC gauge configurations and domain-wall quark propagators used in 
these calculations. The subscript / denotes light quark (up and down), and s denotes the strange quark. The 
superscript dwf denotes the bare-quark mass for the domain-wall fermion propagator calculation. The last 
column is the number of configurations times the number of sources per configuration. For the ensembles 
labeled with "P±A", propagators that were periodic in the temporal direction were computed in addition to 
those with anti-periodic temporal boundary conditions. 



Ensemble 
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bm s 


; dwf 

bm l J 


bm s w f 


# of propagators 


2064f21b676m007m050 


0.007 


0.050 


0.0081 


0.081 


1038 x 24 


2064f21b676m010m050 


0.010 


0.050 


0.0138 


0.081 


768 x 24 


2064f21b679m020m050 


0.020 


0.050 


0.0313 


0.081 


486 x 24 


2064f21b681m030m050 


0.030 


0.050 


0.0478 


0.081 


564 x 20 


2896f2b709m0062m031 


0.0062 


0.031 


0.0080 


0.0423 


1001 x 7 


2896f2b709m0062m031 P±A 


0.0062 


0.031 


0.0080 


0.0423 


1001 x (1+1) 


2864f2b676m010m050 


0.010 


0.050 


0.0138 


0.081 


137 x 8 


2864f2b676m010m050 P±A 


0.010 


0.050 


0.0138 


0.081 


274 x (2+2) 



4. Pion and kaon condensation 



The ground state of the n meson systems that are being studied is a Bose-Einstein condensate 
of fixed z component of isospin (and strangeness in the case of kaons). It is of great interest to 
investigate the properties of such systems. Theoretical efforts have used leading order chiral per- 
turbation theory to investigate the phase diagram at low chemical potential [14] and it is important 
to assess the extent to which these results agree with QCD. Our numerical calculations allow us to 
probe the dependence of the energy on the pion (kaon) density, and thereby extract the chemical 
potential via a finite difference. The results using the coarse MILC lattice are shown for the pion 
and kaon systems in Figs. 2 and 3. Also shown is the prediction from tree-level chiral perturbation 
theory, with which we find surprisingly good agreement. This is encouraging for studies of kaon 
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Figure 1: Three pion (left) and kaon (right) interactions determined from the MILC coarse (blue) and fine 
(magenta) lattices plotted . 
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Figure 2: Dependence of the isospin chemical potential on the isospin density, calculated on the coarse 
MILC ensembles. The curves correspond to the predictions of tree level chiral perturbation theory (dashed) 
[14], the energy shift of Eq. (2.1) (solid) and with the three-body interaction removed (dotted). 



condensation in neutron stars where, typically, tree level chiral perturbation theory interactions are 
assumed amongst kaons and between kaons and baryons 



5. Summary 

Multi-meson systems (and in general multi-hadron systems) have been investigated using lat- 
tice QCD. The calculations presented here provide a first insight into the nature of these systems, 
but much remains to be studied. Recently, we have started to explore the effects of these condensed 
systems on other observables, looking at how the pion condensate screens the potential between a 
static quark-anti-quark pair [16]. 
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Figure 3: Dependence of the strangeness chemical potential on the kaon density. Details are as in Fig. 1. 
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